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We present a discussion of the linear thermoelectric response of an interacting electron gas in 
a quantizing magnetic field. Boundary currents can carry a significant fraction of the net current 
passing through the system. We derive general expressions for the bulk and boundary components of 
the number and energy currents. We show that the local current density may be described in terms 
of "transport" and "internal magnetization" contributions. The latter carry no net current and are 
not observable in standard transport experiments. We show that although Onsager relations cannot 
be applied to the local current, they are valid for the transport currents and hence for the currents 
observed in standard transport experiments. We relate three of the four thermoelectric response 
coefficients of a disorder-free interacting two-dimensional electron gas to equilibrium thermodynamic 
quantities. In particular, we show that the diffusion thermopower is proportional to the entropy per 
particle, and we compare this result with recent experimental observations. 
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I. INTRODUCTION 

It has long been realized that surface (boundary) cur- 
rents provide a significant contribution to the thermo- 
electric response of electronic systems in a quantizing 
magnetic field. The first calculations taking boundary 
currents into account were presented by Obraztsov |^] 
for non- interacting electrons. Such issues have been im- 
portant in recent theories of the thermoelectric properties 
of two-dimensional systems in the integer quantum Hall 
regime. Calculations of the thermoelectric response of 
non-interacting electrons in this regime have been pre- 
sented by a number of authors using various approaches 
Measurements of the intrinsic, "diffusion-and- 
drift thermopower" in the integer quantum Hall regime 
are consistent with theories for non-interacting electrons, 
once disorder is introduced. (For a review of experi- 
ments and theories of the intrinsic and phonon-drag ther- 
mopower for non-interacting electrons, see Ref. 1^.) 

Recently, however, there have been reports of mea- 
surements of the diffusion thermopower in the fractional 
quantum Hall regime It is clear, both from these 

measurements and from what is known of the fractional 
quantum Hall effect, that interactions must play an im- 
portant role in determining the transport properties in 
this regime. 

In this paper we discuss the thermoelectric response of 
an interacting electron gas, paying particular attention to 
the importance of the boundary currents. In section II 
we restate the general expressions for the linear response, 
following an approach first proposed by Luttinger [pd| , 
and since discussed for the quantum Hall regime by Oji 
and Streda [|| . We extend this analysis by deriving gen- 
eral expressions for the local energy-current and number- 
current distributions in gradients of temperature and 
chemical potential. We argue that, even when electron- 
electron interactions are included, both the number and 



energy currents in the bulk of the sample may be sepa- 
rated into "transport" and "internal magnetization" con- 
tributions. The magnetization contribution causes no net 
current to fiow through the sample. However, it can have 
a significant effect on the local current density. We show 
that Onsager relations may still be applied, in a quantiz- 
ing magnetic field, for the transport currents, and hence 
for the net currents through the sample. (The Onsager 
relations cannot, in general, be applied directly to the 
local current densities.) 

In Section III, we consider a two-dimensional electron 
system in the limit of zero impurity scattering, and we 
derive the forms of various transport coefficients in this 
case. In particular, the thermopower coefficient is shown 
to be equal to the entropy per carrier divided by the 
charge of the carrier, a result first derived by Obraztsov 
for non-interacting electrons. In Section IV, we compare 
these results with recent data, at very low temperatures, 
on p-type samples with Landau-level filling fractions near 

V ^ 1/2 and = 3/2 [||. We find that the data at 

V — 1/2 are consistent with an interpretation based on 
a model of spin-polarized "composite fermions" , with a 
reasonable value of the effective mass, but this does not 
seem to be the case at v = 3/2. 

Many of the results of Section II, particularly for the 
net currents, have been obtained previously by Oji and 
Streda |^ , at least for the case where the gradients of the 
potentials and the temperature are constant throughout 
the sample. Many details were omitted from their presen- 
tation, however, and many of the underlying assumptions 
were not stated explicitly. Because there are a number 
of subtle points in the derivation, because there appears 
to have been some confusion in the literature p^ , and 
because the results are of fundamental importance, we 
give here a detailed and general derivation. 
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II. GENERAL EXPRESSIONS FOR LINEAR 
RESPONSE 

A. General Considerations 

We begin by reviewing the "hydrodynamic" assump- 
tions inherent in any theoretical discussion of transport 
coefficients such as the thermal and electrical conduc- 
tivities, thermopower, etc. We restrict our attention to 
small deviations from thermal equilibrium, in samples 
which are very large compared to atomic distances or 
other microscopic length scales, and we shall investigate 
the response to weak external perturbations which vary 
slowly in space and in time. 

We assume here that particles interact only via short- 
range forces, deferring until subsection G the modifica- 
tions necessary in the presence of long-range Coulomb 
interactions. 

The fundamental hydrodynamic assumption is that 
there exists a microscopic relaxation rate Tm, such that 
for perturbations which vary on a time scale slow com- 
pared to Tm , the system relaxes to a state that is close to 
local thermodynamic equilibrium, and where all proper- 
ties of interest may be described in terms of an expansion 
about local equilibrium. More particularly, one identifies 
a set of conserved quantities, which in the systems of in- 
terest to us are the energy E and particle number N, 
and one defines corresponding conserved densities, such 
as the energy density e(r) and particle density n{r). On 
time scales large compared to Tm, one assumes that all 
physical quantities localized near a point r relax to val- 
ues which are determined by the values of the conserved 
densities and their low-order spatial derivatives in the 
vicinity of r. On the other hand, one cannot in general 
assume that the conserved quantities themselves relax to 
their equilibrium values in a microscopic time scale. The 
conservation laws relate the time derivatives of conserved 
quantities to the divergence of associated transport cur- 
rents, and these time derivatives may be very small if the 
length scale of the system is large. 

In systems with short-range forces, the slowest relax- 
ations are typically characterized by a diffusion coeffi- 
cient D, so that the slowest relaxation time for the con- 
served densities is given by tm ~ L'^/D, where L is 
either the size of the system or the wavelength of the 
perturbation, whichever is shorter. Clearly if L is very 
large, tm may be very much larger than t„i. Although 
the overall response to an external perturbation may be 
quite different in the limits where the frequency is large 
or small compared to t^/ , the hydrodynamic equations 
themselves are assumed to apply for time scales r large 
compared to r,„ , regardless of whether r is large or small 
compared to tm- 

Our central focus will be on the particle current den- 
sity J(r) and the energy current density J^{r). At least 
in cases where there are only short-range interactions be- 
tween particles, the hydrodynamic assumptions imply, in 



particular, a locality hypothesis for J and J , viz. that 
J and are determined by the values of e(r) and n{r) 
and their variations, only in the immediate vicinity of 
point r. The currents may also depend on local material 
parameters such as the local chemical composition, im- 
purity concentration, etc., and on the applied magnetic 
field B. We assume the material parameters to be in- 
dependent of time, but they may depend on position in 
cases of interest. If the material parameters are inde- 
pendent of position, then the locality hypothesis implies 
that for variations on a time scale slow compared to the 
microscopic relaxation rate t^^, the current densities J 
and J at point r may be considered to be functions of 
e, rt, and their gradients at point r. 

The locality assumption, central to any hydrodynamic 
description, is difficult or impossible to prove rigorously 
under general conditions. One important piece of the 
physics, which enters the case of quantum systems in a 
strong magnetic field, is that a sample can have non-zero 
"magnetization currents" even in a situation of thermo- 
dynamic equilibrium. In this proof that the elec- 
tron number current J{r) is independent of conditions 
far from r is equivalent to proving that the equilibrium 
magnetization M{r) (defined below) depends only on 
local conditions. For a small sample, M{r) may in fact 
be rather sensitive to conditions at relatively large dis- 
tances from r. For example, for a small metallic loop, at 
very low temperatures, the equilibrium current depends 
on the magnetic flux through the loop, modulo units of 
the flux quantum hc/e, because of the Aharonov-Bohm 
effect. Such effects, however, become negligible in the 
"thermodynamic" limit of large sample sizes. Since the 
equilibrium magnetization density can be related by ther- 
modynamics to a derivative of the free energy with re- 
spect to the applied magnetic field [cf. Eq. (^) below] 
a proof of the locality of M reduces to proving that a 
system has a well-behaved thermodynamic limit for the 
free energy, at any non-zero temperature. 

In our analysis, it will be convenient to eliminate n{r) 
and e(r), in favor of two suitably defined "statistical" 
fields: a local chemical potential /i(r) and a local tem- 
perature T(r). We shall also introduce shortly, external 
"mechanical fields": a potential <j){r) which couples to 
the number density, and a fictitious "gravitational po- 
tential" tp{r) which couples to the energy density. We 
also define an electrochemical potential £^ = jj, + (p. We 
shall see that in thermodynamic equilibrium (even if the 
material parameters vary in space) ^ and (1 -I- "0)2^ are 
constants in space. 

In cases where there is time-reversal symmetry (hence 
B — 0), the currents and J must vanish in thermo- 
dynamic equilibrium. Therefore, the first terms in the 
gradient expansions for and J must be proportional 
to and VT (assuming ■0 = 0). For a quantum me- 
chanical system in the presence of an applied magnetic 
field, however, there may be non-zero circulating cur- 
rents even in a situation of thermodynamic equilibrium. 
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as was noted above. We shall find it convenient to break 
the currents J and into a "transport" part and a 
"magnetization" part according to 



J(r) Jtr(r) + Jmag{r), 

J^{r)^jl{r) + Jl(r), 



(1) 
(2) 



where Jtr and jf^ vanish in thermodynamic equilibrium 
and 



Jmag(r) = V X M"(r), 
jLAr) = V X M^(r). 



(3) 
(4) 



The "magnetization densities" M^{r) and M^{r) are 
defined to be functions of the temperature and chemical 
potential only at the given point r. These functions, in 
turn, may be computed in thermal equilibrium; i.e., we 
may compute the values of {r) and M^{r) assuming 
that fi and T are independent of position. (Any applied 
"mechanical potential" may also be taken independent of 
r.) 

We make a number of observations about the magne- 
tization currents. 

1. If we consider a homogeneous sample with sharp 
boundaries, in thermodynamic equilibrium, then the 
magnetizations and are uniform and the mag- 
netization currents vanish in the interior of the sample. 
However, there will in general be currents flowing on the 
surface of the sample. As one can readily derive from 
Eqs. (^) and (Q), the surface current densities at a point 
T on the boundary are given by 



I = X n, 

= X n. 



(5) 
(6) 



where n is the unit vector normal to the surface, pointing 
outward from the sample. The same expressions for the 
boundary currents are also valid for an inhomogeneous 
sample, in which case the magnetizations , vary 
in the sample. Their values in Eqs. (^) (^) should be de- 
termined at a point inside the boundary close to t. 

We are assuming here, and throughout this paper, 
that the material exterior to the sample is either a vac- 
uum or an "ideal" material with no magnetization of its 
own. Otherwise, we would have a second contribution to 
the edge current from the magnetization of the exterior 
medium. 

2. For a two-dimensional conducting layer in a semi- 
conductor system, the magnetizations M and are 
normal to the layer. If the magnetization lies in the pos- 
itive z direction, then the boundary currents will point 
parallel to the sample edge in the counter-clockwise di- 
rection, looking down at the sample. 

3. The integrated boundary currents given by Eqs. (|^) 
and (|^) depend on the magnetization at a point just in- 
side the sample, but are independent of such details as 
whether the boundary is sharp or diffuse on the atomic 



scale, the concentration of impurities near the boundary, 
etc. In the case of thermodynamic equilibrium, in a uni- 
form sample, where the magnetizations are independent 
of position, the surface currents are divergence-free. This 
is of course necessary since there should be no bulk cur- 
rents in this case. Alternatively, we see that the condition 
of divergence- free surface currents, (|^), (^), for arbitrary 
surface treatments, requires that the magnetizations are 
truly bulk properties, independent of any details of the 
surface. 

4. In our discussion of electron systems, the only par- 
ticles which are allowed to move over macroscopic dis- 
tances are the conduction electrons. The electrical cur- 
rent J*^ is related to the particle current J of the elec- 
trons by J'^ = —eJ, where (— e) is the electron charge. 
For a sample at equilibrium the magnetization M is re- 
lated to the conventional magnetic moment density M, 
in Gaussian units, by 



M ^ {-e/c)M 



N 



(7) 



To avoid confusion, we note, that the quantity M, de- 
fined by Eqs. (|l]) and (|^), has a direct physical meaning 
in terms of the magnetic moment per unit volume only 
if the sample is at equilibrium. In the non-equilibrium 
case, the total magnetic moment is determined by 
the total current density distribution J, including both 
the "transport" component J'tr and the "magnetization" 
component Jmag, as given by the general formula 



M 



d^r[rx (J(r)- J)], 



(8) 



where the overline denotes the volume average. More- 
over, since Jtr cannot, in general, be expressed as the 
curl of a vector field, the magnetic moment of the sam- 
ple in the non-equilibrium case cannot be expressed as 
an integral of any local magnetization density. 

5. For a uniform macroscopic sample, in thermal equi- 
librium, the magnetization may be related to other ther- 
modynamic quantities. By definition, the magnetization 
M is equal to —{l/V)dE/dB, when the entropy S, the 
volume V, and the electron number N are held fixed. 
More generally, we may write 



TdS ^dE + PdV + MV ■ dB ~ ^idN, 



(9) 



where /i is the chemical potential of the electrons, and P 
is the electron "pressure" . (We are assuming here short- 
range forces between the electrons.) From the extensivity 
properties of S*, E and iV, it then follows that 



n/i = e — Ts + P, 



(10) 



where e and s are the energy and entropy per unit vol- 
ume. Then using (|^) we find 



N 



M 

P = 



c dP 
e dB 



/ n{^,\ T, B) dfi . 

J — oo 



(11) 

(12) 
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6. If there are temperature or electrochemical poten- 
tial gradients, there may be non-vanishing magnetiza- 
tion currents in the interior of an otherwise uniform sam- 
ple. Alternatively, if there are gradients in the material 
parameters, bulk magnetization currents can be present 
even at thermodynamic equilibrium. 

7. Under all conditions, the magnetization currents 
J mag and J^ag divergence-free. As a consequence, 
these currents do not make any contribution to the net 
current flows that are measured by conventional trans- 
port experiments. Specifically, consider any closed curve 
C that encircles the sample but is exterior to it, and let 
S he a, surface spanning this contour. The total magne- 
tization currents crossing the surface must be zero, by 
Stokes' theorem, and the total currents Is and Ig cross- 
ing the surface are obtained by considering the transport 
currents alone: 

Is= [ h-J(fS= [ h Jtr (fS, (13) 
J s J s 

/|= ( fi-J^ <fS= [ h-Jtr^d'^S, (14) 
Js Js 

where h is the local normal to S. In a d.c. transport 
experiment, where V • J = V • = 0, the currents 
Is and Ig will be independent of the particular surface 
chosen. This argument applies equally well to a singly- 
or multiply-connected sample including the case where 
C threads a hole in the sample, as in Figure 1. In the 
case of a two-dimensional sample, the surface S becomes 
a curve traversing the sample, while Is and Ig are the 
total currents across the curve. 

Current Current 
Lead Lead 




FIG. 1. Contour C, external to sample, is spanned by sur- 
face S. The total current through 5* is equal to the transport 
current across 5", as the magnetization current gives no contri- 
bution. The illustration shows a multiply-connected sample. 

In the remainder of this section, we shall apply the 
above considerations to the general hydrodynamic de- 
scription of an electron system in a strong magnetic field. 
Our goal is to obtain some relations among the various 
transport coefficients - both for the transport currents 
Jti{r) and jf^{r) and for the total local currents, J(r) 
and J^{r). We also relate the transport coefficients to 
microscopic expressions involving correlation functions 
for currents in the equilibrium state. 



Our strategy can be best illustrated by considering a 
simple example where the electrons are subject to a weak 
external potential (t>{r), while the temperature is main- 
tained constant, say by contact with a substrate. Then 
inside the sample we may write 

Jt,(r) = V/i(r-) - L^'^^V<l>{r), (15) 

where fi{r) is the local chemical potential for the elec- 
trons, and N^^') and L^^^ are second-rank tensors, which 
we denote "transport coefficients." In thermodynamic 
equilibrium, the electrochemical potential £^{r) — nir) + 
4>{r) is independent of position, and Jtr, by definition 
is zero. Since this must be true for arbitrary (t>(r), we 
immediately obtain the "Einstein relation" N^^^ = 

Macroscopic equations for the response of the system 
to a time-dependent perturbation <j)[r,t) are obtained 
by combining ( p^ ) with the conservation law {dn/dt) — 
— V • Jtr- (Recall that V • Jmag = 0.) If there is a pe- 
riodic disturbance in ^(r) with a wave vector q\\x, then 
the electron density will relax towards the equilibrium 
state, with ^(r) — const, at a rate 7^ = Dq^ ^ where the 
diffusion constant D is given by 

D = L^:^}/{dti/dn). (16) 

(Again, we assume short-range interactions between the 
electrons, so dpi/dn is finite for g — s- 0.) A microscopic 
expression for L(i) can be obtained by using quantum 
mechanics to calculate the response of the system to an 
infinitesimal time-dependent perturbation (p, applied at 
a frequency uj which is small compared to microscopic 
frequencies t^^, but high compared to 7,, so that the 
density n does not have time to change significantly. One 
thus obtains an expression for L^^^ in terms of a two-time 
correlation function for fluctuations in J in the equilib- 
rium state. 

In order to generalize this procedure to the case of non- 
uniform temperature, we follow the work of Luttinger 
|]lT| which involves the introduction of a fictitious "grav- 
itational potential" il){r), which enters the Hamiltonian 
through its coupling to the energy density. If linear re- 
sponse to ■(/; is calculated, the response to a temperature 
gradient may be obtained from the Einstein relations. 
Our situation is more complicated than the case consid- 
ered by Luttinger, however, because the Einstein rela- 
tions apply only to the transport currents, and not to 
the total currents. 



B. Current Operators in the Presence of Electrical 
and Gravitational Fields 

The linear response to mechanical fields (electrostatic 
and gravitational) in a quantizing magnetic field has been 
presented previously H. However, since it is important 
to the rest of our discussion, we shall review this here. 
We consider a Hamiltonian of the form 
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N 



N 



^[p^ 



■ eA{r,)/c\ 



2m 



where A{r) is the magnetic vector potential, V{r) is the 
scalar potential energy including the confinement and 
disorder potentials and the periodic potential of atomic 
cores, and u.^j describes the interparticle interactions. 
Following Luttinger, we introduce the number and en- 
ergy density operators ||ll| , [l3|| 



i 



(18) 
(19) 



where Si = S{r — ri), hi is defined in equation (|17[), 
and the curly brackets indicate the anticommutator, 
{A,B} = AB + BA. 

We are interested in the response of the current-density 
to time varying external "electrostatic" and "gravita- 
tional" potentials (/) and -01 respectively. These fields 
couple to the number and energy densities according to 
the Hamiltonian 



Ht = J d-'rhrir), (20) 
hrir) = h{r) + (/)(r)p(r-) + ii;{r)h{r). (21) 



[We call the function ip{r) defined by Eqs. (|20|) and (|21 
the "gravitational potential" to follow the terminology 
of the original work by Luttinger ||ll| , although the true 
gravitational potential would also be coupled to the mass 
density rather than just to the first-order term of the rel- 
ativistic energy expansion.] 

The conservation laws for energy and particle number 
imply that the Heisenberg equations of motion for p and 
hx, under the Hamiltonian Ht, may be written in the 
form 



J 




(17) 



r 



dt 
dhxir) 
di 



--[p,HT]=-V-J{r), 



~dt 



-V-J^{r) 



dhrir) 
dt - 



(22) 
(23) 



where J and are here operators for the particle and 
energy currents. The second term on the right-hand side 
of ( P3[) occurs when there is a time dependence to (/)(r) 
or i^{r). The above equations constrain, but do not de- 
termine uniquely, the forms of the operators J and J , 
as we may in principle add to them an arbitrary function 
whose divergence is identically zero. A requirement of 
strict locality for the particle current, however, together 
with Eqs. ( po| ) and (pi|), imposes the form 



J{r)^j{r)[l + ^j{r% 
where j{r) is the particle current for ijj ~ 0: 



3A{ri)/c]/m. 



(24) 

(25) 
(26) 



The requirement of strict locality cannot be applied ei- 
ther to the energy current or the energy density in a non- 
relativistic theory with interactions of non-zero range. 
However, the form of the energy current is adequately 
restricted for our purposes if we require that it depends 
only on the positions and velocities of particles in a small 
neighborhood of r. Although various forms of the en- 
ergy current may still be written down, which are consis- 
tent with this requirement and with a requirement that 
Eq. ( p3| ) be satisfied exactly on the microscopic scale, we 
shall adopt here the particular definitions 



J^(r) = j^(r) + 4'{r)j{r) + 2VXr)j^(r), 

(r) ^ i ^{/i,, {vr,S,}} + i ^ ^{(«7 + v]), r^^{r)}, 

i i^j 7 

r;^^(r)^(rf-r^")^;] A,,(r), 

^ii(^) = j ds 5[r - Tj - {ri - rj)s], 
Jo 



(27) 
(28) 

(29) 
(30) 



where = —duij/drf is the force on particle i due 
to its interaction with j. We note that by definition 
is the energy current associated with the unperturbed 
Hamiltonian Hq, and equation (|2^) is valid to first order 
in -0, with the requirement that <f) and ijj vary very slowly 
in space compared to the range of the interaction Uij . It 
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is not difficult to show, under these circumstances, that 
exactly satisfies the Heisenberg equation of motion 
(p3|). To demonstrate this, it is convenient to integrate 
Eq. (^) over a small volume 5V , and to write the inte- 
gral of V • as an integral over the surface S of the 
volume. We also note that 



5 



/ (r) {r, -rj)-dS = Xj - Xt 
Js 



(31) 



where Xi = if particle i is inside the volume enclosed 
by 5, and zero otherwise. 

With our definitions, the contribution to the energy 
current from the interaction potential Uij is distributed 
along the straight-line segment joining and rj, sim- 
ilar to the distribution of the ener gy current adopted 
in Ref. By contrast, in Ref. Luttinger em- 

ployed a different form for , where the energy current 
is supposed to be localized at the two points and rj. 
Luttinger's expression is only approximate on the micro- 
scopic level, since it does not strictly satisfy Eq. (23). 
Our more precise formula and Luttinger's are equivalent, 
however, when integrated over any volume large com- 
pared to the range of Uy . In fact, in the present paper 
we have no real need for a precise form for j^{r), but 
will only make essential use of (E^. 



C. Thermodynamic Equilibrium 

We first consider the situation of thermodynamic equi- 
librium. A grand-canonical ensemble, for the perturbed 
Hamiltonian Ht may be described by a density matrix 



w 



(32) 



where ^ and /3 = T^^ are Lagrange multipliers, which we 
may describe, respectively, as the "electrochemical po- 
tential" and the "inverse thermodynamic temperature." 

The properties of w are particularly simple in the case 
where ip and (j) ctre constants, independent of position. 
Then the eigenstates of Ht are the same as the eigen- 
states of Hq, and the density matrix w is identical to a 
density matrix wq , of the form 



Wo 



with 



T-i=/3(l + V'), 

/i = (e-<^)/(i + v)- 



(33) 



(34) 
(35) 



We denote T as the "internal temperature," and /i as the 
"internal chemical potential" . The internal energy den- 
sity e(r) = {h{r)) and the particle density n{r) = {p{r)) 
are the same in the two cases, as is the entropy S = 
—tT(whiw). Note that e and n may depend on position, 
even in equilibrium, if the materials parameters vary in 
space. 

The particle and energy currents are not the same in 
the two cases, however, because the operators themselves 
are modified by the presence of tp and (j). Using ( |2^ ) and 
(p7|), we see that 



{J{r)) = {j{r))+^j{j{r)), 



(36) 
(37) 



The expectation values on the right-hand side represent 
the particle and energy currents in the "unperturbed" 
state, with (j) — ip = 0, temperature T, and chemical 
potential /i. [Note: in Eq. (|3^), as in Eq. (|2^ ) we have 
dropped terms of order ip'^, since we always consider ijj 
to be small.] 

The currents in Eqs. (^ ) and ( p7| ) are purely mag- 
netization currents, since the transport currents vanish, 
by definition, in equilibrium. In the case of a uniform 
sample, the magnetization currents are confined to the 
boundaries. More generally, they are given by Eqs. (H) 
and (^). This implies that in the equilibrium state, with 
constant (j) and ^Z", we have, to first order in the pertur- 
bations, at any point r in the sample 



N/ 



(38) 
(39) 



where Mq and Mq are the unperturbed magnetiza- 
tions, (corresponding to '0 — 4' ~ 0); which we have 
written explicitly as functions of the chemical potential 
/i and internal temperature T. (The functions and 
Mq are also implicitly functions of the material parame- 
ters, such as chemical composition and impurity concen- 
tration, in the neighborhood of the point r.) 

The results (^,^,^,^ can be generalized to the case 
when (j) and ip vary in space. Equation (|3^) is still valid 
in this case. If the length scale of the variation is suffi- 
ciently large, the entire system can be considered, within 
the standard approach, as consisting of small subsystems 
weakly interacting with each other at imaginary borders. 
Each of these subsystems, if allowed, will reach a local 
equilibrium state which depends on the local values n{r) 
and e(r). We assume, of course, that the magnetic field 
is constant. One can then introduce local thermody- 
namic parameters (chemical potential, temperature, en- 
tropy density, etc.) which are functions of the two inde- 
pendent variables, n and e. The values of extensive ther- 
modynamic parameters will be then given by the integral 
of the corresponding local densities over the sample area. 
In particular, for the total energy, total entropy and total 
number of particles we have 

Et = {Ht) = j [e(r) + n{r)(j){r) + e(r)V'(r)] dV, (40) 
S = j s(e(r),n(r)) dV, (41) 
N = j n{r) d^r, (42) 

where s(e, n), the local entropy density, is the same func- 
tion of e and n as for a uniform system. The accuracy 
of this approximation is limited by terms (Ve)^, (Vn)^ 
which reflect the interaction between adjacent subsys- 
tems. 

It is more convenient to choose, instead of e and n, a 
local "internal" temperature T{r) and chemical potential 
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/i(r) as independent thermodynamic variables. Following 
Luttinger, we define these parameters in the same way 
as it is conventionally done for a homogeneous system 



de 



dn 



(43) 



[Of course, these definitions are consistent with the ther- 
modynamic relation (H), interpreted in a local sense, with 
V and B fixed.] 

The equilibrium state defined by the density matrix 
(^) is the state which maximizes the functional $ = 
S - (3{Et ~ ^N). Using @, (|l|) and and setting 
the variational derivatives of equal to zero, we find 



Kr) = [e-0(r-)]/[l + ^(r)] 



(44) 
(45) 



which is the generalization of (^^ and ( p5| ) to the non- 
uniform case. The equilibrium conditions may be alter- 
natively written as 



5(m/t) 

5{l/T) : 



-4>{r)/To, 



(46) 
(47) 



where Tq"^ — (3 is the unperturbed temperature of the 
sample. If To is replaced by the perturbed temperature 
T, these conditions reduce to the conditions for thermal 
equilibrium obtained by Luttinger |pT[| , which are valid 
for small changes in (j) and ip. 

We now consider the magnetization currents. By def- 
inition, the transport currents vanish in an equilibrium 
state, so the magnetization currents are the total local 
currents in that case. If (/)(r), ip{r), T{r) and /x(r) vary 
sufficiently slowly in space, it is clear that {r) and 
M^(r) wiU be still given by Eqs. ( p8| ) and (|39|), provided 
that the arguments /i and T of the functions Mq and 
Mq are evaluated at the position r. (Corrections will be 
at least second order in the gradients, if one is away from 
the boundary of a sample.) Then using Eqs. (||) and (|^), 
for a sample with uniform chemical composition, far from 
the boundaries, we find 



rbulk 

mag 



r-E,bulk 



dMl_ 



dT 
dT 



N 



VT-M^ X Vij), 



X VT 



X 



2M^ X Vi), 



(48) 
(49) 



to leading order in </> and ij}. There will also be addi- 
tional magnetization currents at the boundaries, given 
by Eqs. (I, (|), H) and (H). 



D. Non-Equilibrium States 

We now consider the non-equilibrium situation, where 
Eqs. (46) and ( ^7| ) are not satisfied and where cj), -0, 
T, and /i may be in general time-dependent. The local 
temperatures T{r) and /i(r) are defined to be the same 
functions of €{r) and n{r) as in the equilibrium case. We 
continue to define the magnetization currents and mag- 
netizations by (38), (^, ( p8| ) and (|9|), and we define the 
remaining contributions to the currents to be the trans- 
port contributions. 

By the locality hypothesis, the transport currents must 
be given by a sum of terms proportional to the gradients 
V(^, VV', VT and V/i, in the limit where the applied 
fields are small. Our aim is to determine the coeffi- 
cients of these terms. The requirement that the trans- 
port currents vanish when (^6|) and (^^ are satisfied 
means that the gradients enter only in the combinations 
+ TV(/i/T)] and [V^ - TV(1/T)]. This observa- 
tion is the generalization of the "Einstein relation" to the 
case in which there is a non-uniform t/;, and allows the 
response to the statistical fields /i and T to be related to 
the response to the mechanical fields (j) and V'- 



To determine the remaining unknown coefficients, we 
consider a particular situation, where the potentials 4> 
and vary periodically in time, with a characteristic fre- 
quency bj small compared to the microscopic relaxation 
rate but large compared to the size-dependent relax- 
ation rates r'j^^ for the energy and particle density. (This 
is the "rapid case" of Ref. |ll).) In this situation the lo- 
cal values of e and n are not changed from their initial 
values, so that and T remain constant throughout the 
sample. 

We concentrate on a region far from the boundaries. 
Since VT = = 0, we may write the total currents as 

jbuik ^ 2(i)(_V(/.) + L(2)(-Vi^), 



rB.bulk 



i(3)(„V(/))-f l(4)(-vv;), 



(50) 
(51) 



with transport coefficients that can be expressed in terms 
of time-dependent correlation functions in the framework 
of the standard Kubo formula [pi]] : 

1 foo j-fJ 
= 1™, V j,'^^ jf ^"^ ~ ' ^^^^ 

L(2) ^ lim - / dt e-* /d/3'(jj^(-t - z/3')jo"(0)), (53) 

ii^) ^ lini _ Jdt e-^* jd(3'{foi-t - ^/3')jo ''(O)), (54) 

1 roc i-fS 
LW = lim- /die-^* /d/3'(jr(-t-z/?')jo"^(0)). (55) 
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Here the subscript on each time-dependent current op- 
erator Jq and Jq indicates the q —> Umit of its spatial 
Fourier transform, and the angular brackets denote the 
quantum mechanical and thermodynamic average in the 
equilibrium state of an infinite system at temperature 
(3^^ = T and chemical potential fi. 

In order to find the transport currents when VT = 
W fj, — 0, we subtract from (pOh and (|5l|) the bulk magne- 
tization currents, given by (Eq) and ([49|). We thus find 

Jtr = -L[^W(I) - [i^^^ - • e] VV', (56) 
Jf, = _ ■ e)V0 - - 2M^ ■ e] (57) 

where e is the unit antisymmetric three-tensor: (e)a/37 = 
1(— I) if the a(3"f is an even (odd) permutation of xyz, 



and zero if two indices are equal. By definition, transport 
currents flow only in the bulk of the sample. To simplify 
notation, we have therefore dropped the label "bulk" in 
the above equations, and will do so for all subsequent 
transport current densities. 

Now we can consider the general case, for an arbitrary 
value of lutm, where V/i, VT, VtA, and V0 may be 
all non-zero. As remarked earlier, the condition that the 
transport currents vanish in equilibrium requires that the 
gradients enter only in the combinations [V0-|-TV(/x/T)] 
and [VV' — TV(l/r)]. (In other words, the Einstein re- 
lations, extended to the case where VT and are dif- 
ferent from zero, apply to the transport currents in the 
presence of finite magnetic field.) Thus we have, in the 
general case: 



I 

Jtr = -L(i)[V0 + rV(/i/T)] + (L(2) _ ■ e)[- + TV(l/r)], (58) 
Jf^ = - ■ e)[V(j) + TV(^/T)] + {L^^'> - 2M^ ■ e)[-V^j + TV{1/T)]. (59) 



Of course, (|58| ) and ( p9| ) reduce to the Luttinger for- 
mulae 11 1 when the magnetic field is absent, so that 
= M^ = 0. 

To obtain the local currents, in a uniform sample far 
from the boundaries, in the general case where VT, V/i, 
V0 and V?/' are all independent, we must add the magne- 
tization currents, given by ( ^8|) and ( |49| ) to the transport 
currents given by ( ^8|) and (|59|). In addition, there will 
be boundary currents from the magnetization. In the ab- 
sence of the fictitious gravitational potential (tp = 0) the 
change in the integrated boundary currents at a point r, 
Eqs. (^) (^, is determined by the changes in the tem- 
perature and chemical potential at that point, according 
to 



61 : 



= -n X 



QMS 
dMj 



N 



■5ii 



dT 



N 



■6T 



dT 



N 



(60) 



(61) 



Although we have introduced the fictitious field ip to de- 
rive linear response, we will have no further need of it. 



and hence shall set ip = for the remainder of this paper. 

It is convenient to introduce the electrochemical po- 
tential, defined by 



ar) = ^l{r) + (f>ir). 



(62) 



An ideal voltmeter, (with leads that have no ther- 
mopower) will measure the difference in ^ between two 
contact points. In a thermodynamic equilibrium state, 
the value of ^, as well as the temperature T, will be con- 
stant throughout the system. 

It is also convenient to define a transport heat current 
density as 



tr — tr tr- 



Then we may rewrite ( p8D and (|59|) in the form 
= -N^^^Vi - N^'^\VT)/T, 

where 



(63) 



(64) 
(65) 



/^(4) = ^(4) 



m(x(3) 



M 
M 





1(2)- 



e. 



2(M^ 



(66) 
(67) 
(68) 
(69) 



Although our derivation has assumed (p to be infinites- 
imal, the final results for the transport currents, given 
by (|6^)-(|69[) and (52)-(|5^) are written in a form that 
remains valid when (f> is not small. This accounts for 
the slight differences (which are beyond leading order in 



the driving fields) between the above coefficients, N^^\ 
and those obtained by direct substitution of Eqs. ( [S^ ) 
and (^9|) into Eq. (p3|) . When (j) is not small, the ex- 
pressions ( ^0| ) and ( |51| ) should contain additional higher- 
order terms in 0; the coefficients L appearing above can- 
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not therefore be defined by these expressions, but are as- 
sumed to be defined by the Kubo formulae (|5^,^,^,^ . 
In the case of finite equation ( ^9| ) for the bulk mag- 
netization current J^'il"''^ must also be modified by the 



mag 



addition to the right-hand side of a term J^ag i the edge 
current is still given by and (^). 



E. Onsager Relations 



J?,, given by 



The transport coefficients for Jtr and 

(|6^)-(|69|), obey Onsager symmetry relations of 
the form 

iVW(B) = ivW(-B), (70) 

Ni^^{B) = N('2{-B), (71) 

n!^]{B) = n(^2{-B). (72) 
To see that this is the case, we first establish that the co- 



efficients L 
same form 

relating La\ 



'ay ^ — 1, 2 , 

as 



3, 4) obey Onsager relations of the 
. This follows from the expressions 
to the current correlation functions, (|5^)- 
(^5|), and the invariance of the Hamiltonian Hq under 
simultaneous reversal of time and magnetic field. Sec- 
ondly, we note that the magnetizations Ai'^ and 
reverse sign under reversal of B. Equations ([70|)-(|7^) 
follow directly. 

In contrast the local current densities do not satisfy 
the Onsager relations. The local current densities dif- 
fer from the local transport current densities by the bulk 
"magnetization currents" . These additional magnetiza- 
tion contributions to the local response in general depend 
on which driving field is applied (e.g. whether it is V/i 
or V0), and give rise to a set of transport coefficients, 
one for each driving field, which differ by terms propor- 
tional to the gradients of magnetization densities. Since 
no general symmetries relate these additional magnetiza- 
tion terms, neither the Einstein relations nor the Onsager 
relations hold locally for the total current response. 

As a specific example, consider the local bulk currents 
J and J'^ = — in terms of V/i and VT, under 
conditions where W(t> — Wtp = 0. By combining (|6^)- 
(|^), with (0), (|), (H) and (||), one may readily obtain 
expressions for the appropriate coefficients Ca 



£(2) ^ ^(2) ^j. 



£(3) ^ ^(3) ^ 



N 



dT 



N 



(73) 
(74) 



As far as we are aware, there is no symmetry relating the 
derivatives of Mq and with respect to and T, 

so these coefficients apparently do not satisfy the usual 
Onsager relations. 

Moreover, if we consider instead a situation where 
V/i = VT = 0, but Vcj) ^ 0, we find according to (|5C 



and (|T]), = -(L(3) _^l(i))V0. This coefficient, 
which is different from C'^'^^ , is also clearly not related to 
by Onsager symmetry in the general case. 

Although the contribution of the magnetization cur- 
rent does not appear in standard transport experiments, 
the local electric current distribution is measurable, at 
least in principle, by a sufficiently sensitive measurement 
of the magnetic field generated by currents in the sample. 
Since the magnetization current does not dissipate Joule 
heat, it could not be detected using the well-known tech- 
nique based on local luminescence intensity ijl^. Mea- 
surements of the local electric field distribution by means 
of the electro-optic effect @ do not serve this purpose 
either, since the magnetization current is related to the 
"statistical fields" V/j, and VT rather than to the electric 
field (l/e)V0. In the case of the energy current or heat 
current, we are not aware of any reasonable method for 
direct measurement of the local currents. 

We illustrate our results further with the example of 
the number current J for a two-dimensional electron sys- 
tem with no disorder, at a Landau-level filling fraction v* 
for which the electron system is incompressible at T 0. 
At the filling fraction v* , we have n — v*e\B\/hc. Then 
we have 



dM 



N 



c dn 
"e dB 



-—z sgn(B,), (75) 



where the last equality follows from the incompressibil- 
ity condition for T — > 0. Here 1; is a unit vector directed 
upward from the plane. 

The coefficient i*^^^ for this system is given by the 
Hah conductivity: i'^) = (^^^ sgii{B,^){z ■ e). On the 

other hand, substituting £^ — + cj) into Eqs. ([T]), ( |4^ ) 
and (|5^), we see that the coefficient describing the lo- 
cal response to V/i vanishes, and all number current 
fiow in the bulk is due to the electric field: J^^^^ — 
/h)sgn(Bz)[z X V0]. The current driven by an in- 
homogeneous chemical potential is localized at the edge, 
and is given by ^7 = {i/* /h)sgn{Bz)6fi{r)[h x z], where 
n is a unit vector in the plane, perpendicular to edge, in 
the outward direction. 



F. Inhomogeneous Samples 

Although we have concentrated so far on the case of 
a homogeneous sample with boundaries that are sharp 
compared with the overall length scale, it is easy to gen- 
eralize our results to the case of a sample whose mate- 
rial parameters, such as chemical composition, vary on 
a macroscopic length scale. The formulae (^) and ( |65| ) 
still hold locally for the transport currents in this case, 
with the qualification that the transport coefhcients A^^*-* 
depend on the local material parameters, and can there- 
fore vary from one place to another in the sample. The 
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magnetization currents at the boundary of the sample are 
still given by ( |60| ) and (|6l| ) , with the qualification that the 
functions Af g and Mq may also vary from one place to 
another because of their implicit dependence on the local 
material parameters. The magnetization currents in the 



bulk of the sample are no longer given simply by Eqs. ( |4q ) 
and (^) however. If we denote the material parameters 
by a set of variables {rji}, and we set Vip = 0, then ( |4^ ) 
and (^) should be replaced by 



N 



rbulk 



rB.bulk 



9Mg 

dMl_ 



X Vfj. 



dm 



X V?7j 



X v.; 



bulk 
mag • 



(76) 
(77) 



[These expressions follow directly from (||) , fl) , ( |38| ) and 
(^).] In the quantities OMq /djji and c^Afg /i9?7i, which 
appear in (^6|) and (|77|), it is necessary to keep not only 
the zeroth order terms, but also the first-order changes 
engendered by the deviations 5^{r) and 6T{r). 



G. Long-Range Forces 

The derivation in the previous subsections, and vari- 
ous intermediate results, require modification when there 
are long-range forces due to unscreened Coulomb inter- 
actions. In two-dimensional structures, this is the case 
when the metallic gate is either absent or is situated far- 
ther than the characteristic length scale of the fields ap- 
plied. 

If uncompensated electric charges are present, there 
can be energy transport over long distances via the 
macroscopic electric field, and the total energy current 
at a given point of space does not, in general, depend 
solely on the state of the particle system in the immediate 
neighborhood of that point. Moreover, the convention of 
Eq. (p8|), in which the interaction contribution to the 
energy current is concentrated along the line-segments 
joining each pair of particles, is not generally used in this 
case. 

The most convenient approach is to break up the in- 
terparticle interactions into a long-range piece, mediated 
by the macroscopic electric field E{r), and a short-range 
piece, which includes everything else. The macroscopic 
field is supposed to be averaged over a region sufficiently 
large that fluctuations in the field, arising from thermal 
or quantum fluctuations in the microscopic electronic 
charge density, can be neglected. Thus there is no en- 
tropy transport via the macroscopic field. As a result, we 
find that with appropriate definitions we can write local 
hydrodynamic equations for the heat current and parti- 
cle current which are similar in form to the equations 
derived for short-range forces. (We assume here that we 
are working at a temperature sufhciently low that heat 
transport via the radiation field may be completely ne- 
glected. Formally, this assumption is imposed by taking 
the limit where the speed of light c is infinite.) 



To make these points clearer, we note that since charge 
is locally conserved, charge fluctuations with the longest 
range effects are electric dipole fluctuations. In the 
absence of the radiation field, the interaction between 
dipoles at two different points falls off as the inverse cube 
of the separation, and the rate of energy transfer due 
to random thermal motion would be expected to fall off 
as the inverse sixth power of the separation. At large 
distances this process will be much slower than the con- 
ventional process of heat conduction (already taken into 
account in our discussions) wherein energy is transported 
diffusively via a series of many short jumps. 

To proceed formally, we redefine the interaction Uij 
appearing in subsection B to include only the short- 
range part of the Coulomb interaction, after effects of the 
macroscopic electric field E are subtracted out. First, we 
choose a truncation radius rcut which is much larger than 
the average inter-particle distance but much smaller than 
the macroscopic scales of external fields and, if the sample 
is not uniform, the scale of the equilibrium density varia- 
tion. We split the Coulomb interaction potential into the 
sum of two terms, u-°*(r) = uf-°'^^{r) + ii[™^(r), where 



.short 



(r) decays rapidly at r rcut, and u-- ®(r) con- 
tains the long-range tail of the interaction and changes 
smoothly at distances r <C rcut- Then we replace by 

^, short ' ' 



in the definition of Hq and hi, Eq. (|17|), and include 

the long-range component u[™^(r) in the self-consistent 
macroscopic field E. 

We continue to define h{r) by Eq. (|l9|), and define 
the internal energy density e{r) as the expectation value 
{h{r)). Then e may be interpreted as the matter con- 
tribution to the energy density. The total energy (with 
■0 = = 0) is then given by 



E= (Pr 



Stt 



\B\ 



Stt 



(78) 



where k is the dielectric constant of the background mate- 
rial, and we assume the background magnetic permeabil- 
ity is unity. The energy current will similarly be broken 
up into two parts 



■ E _ -E 
Jtot ^ J 



cEx B 

47r ' 



(79) 
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where the second term is the standard contribution from 
the macroscopic electromagnetic fields, and , which 
we may think of as the matter contribution to the energy 
current, is defined by Eqs. (|2^)-(|30|), with Uij replaced 
by in the definition of F^". Of course, E and B are 

determined self-consistently from the macroscopic cur- 
rent and charge distributions using Maxwell's equations 
in the static limit. Equations (^8|) and (|79| ) are asymp- 
totically correct in the limit where Tcut is large compared 
to the microscopic scale but small compared to the scale 
of variation of E. 

For a two-dimensional electron system in an external 
magnetic field, the magnetic fields arising from currents 
in the sample are generally very small, and may be omit- 
ted from the term \B\'^/8tt in Eq. Thus this term 
is independent of the state of the electron system and 
may be ignored if desired. On the other hand, the mag- 
netic fields generated by the currents in the sample must 
be included in the second term on the right-hand side 
of (v%, because the speed of light appears as a prefac- 
tor. Note also that e and j are restricted to the two- 
dimensional layer, but the electromagnetic contributions 
to and d?^) extend into the space outside. Below we 
focus on the matter part of the energy current. 

It is now possible to redo the arguments of the previous 
section with little modification. We restrict our attention 
to the situation where the applied magnetic field is inde- 
pendent of time and assume that the macroscopic electric 
field may be derived from a scalar potential 

E{r) = -V$(r). (80) 

The potential $ is obtained self-consistently, and includes 
effects of any macroscopic time-dependent variations in 
n(r), as well as the effects of any static charges present in 
equilibrium and the external perturbations embodied in 
(j). In general, we cannot consider that $ is infinitesimal, 
even in equilibrium. 

In order to repeat our previous derivations, one needs 
to consider a non-zero gravitational potential ip; we 
present here the final result, and hence set ^ = 0. We 
redefine the electrochemical potential as 

C(r) = MW-e<f(r). (81) 

The definition ( p7|) of the energy current in the perturbed 
system is now replaced by 

jE ^ -E _ (-32) 

where J — j is the particle current density, and was 
defined above. As we did in previous sections, we split the 
local currents J and into magnetization and trans- 
port parts, and we define a transport heat current as 
given by Eq. ( |6^ ) . Then J^. and Jtr obey equations iden- 
tical to (||)-(|9|), with L'^al defined by Eqs. (||)-(||) in 
terms of the correlators for the matter currents j{r) and 
j (r) under the Hamiltonian Hq for a uniform system in 



equilibrium, with (j) = ip = 0. The transport coefficients 
Na-^ obey the same Onsager relations as before. 

Note that different definitions of the macroscopic elec- 
tric field, as may be obtained by different choices of the 
truncation radius rcut, will generally cause an exchange 
of contributions between the chemical potential fi{r) and 
the electrostatic potential <I>(r). This will also trans- 
fer contributions between the first and second terms in 
Eq. (|8^), leaving the sum unchanged. As long as 
the different values of rcut remain sufficiently large, the 
change in the truncation should not affect the coefficients 
i'^'-' since the correlators in Eqs. (|5^-(p5|) are sensitive 
only to short-range properties of the system. 

We also note that in realistic two-dimensional systems, 
the component of the electric field perpendicular to the 
electron layer, arising from charges on gates or from ion- 
ized impurities displaced from the layer, may play an 
important role in confining the electrons to the layer. It 
is therefore important to include this part of the macro- 
scopic electric field in the unperturbed Hamiltonian Hq 
and in the definition of the energy current when calcu- 
lating the correlation functions that appear in (p2[)"(p5|). 
Formally this can be done by including the perpendic- 
ular confining field in the one-body potential V which 
enters Eq. (O), and excluding it from the macroscopic 
potential <I>(r). Then the equilibrium $(r) is a constant 
in the direction perpendicular as well as parallel to the 
layer, and can be safely omitted from Hq. In fact, it may 
be convenient to include the entire equilibrium value of 
the electrostatic potential in V, even in an inhomoge- 
neous system, so that ^ describes only long-wavelength 
fluctuations about equilibrium. 



H. Additional Remarks 

Throughout this section, we have treated the fields /z, 
T, ■0 and (p (or $, in subsection G) as independent vari- 
ables, which can be arbitrary functions of position, sub- 
ject to the constraints that gradients are small, and that 
4> and V' E^re infinitesimal. (The fictitious field ip was 
set equal to zero in the latter part of the section.) The 
transport currents depend on gradients of ^ and (j> (or 
<&) only through the combination V^. In a d.c. trans- 
port experiment, there are strong additional constraints 
arising from the requirements that the currents must be 
divergence-free in the interior of the sample, and must 
satisfy appropriate constraints at the boundaries. Typ- 
ically, these conditions completely determine the spatial 
variations of ^ and T throughout the sample interior, 
when boundary values of the fields are specified, or when 
current flows through the boundaries are given. 

For the case of a two-dimensional electron system on 
a three-dimensional substrate, the situation is slightly 
more complicated. We shall be concerned with situations 
where the substrate is an electrical insulator, so that the 
divergence of the two-dimensional electric current is re- 
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quired to be zero in the analysis of experiments. On the 
other hand, we consider the thermal coupling to the sub- 
strate, via absorption or emission of phonons, to be small 
but not zero. Then on length scales large compared to 
{DxTepY^^ , where r^p is the electron-phonon relaxation 
time and Dt an appropriate thermal diffusion constant 
for the isolated electron system, the divergence of J'^ is 
not necessarily zero. Instead, one should take the value 
of T{r) to be an independent variable determined by con- 
ditions in the substrate. 

Although our previous discussions assumed the sam- 
ple to be isolated from its environment except at its 
edges, the transport equations derived above should re- 
main valid provided the electron-phonon coupling is suf- 
ficiently weak that Tep is large compared to the micro- 
scopic times necessary to establish local equilibrium in 
the electron system. 

In order to calculate the transport currents, it is not 
generally necessary to find the separate portions of 
arising from V/i and from the electric field. This is nec- 
essary, however, if one wishes to obtain the local current 
distribution. In practical situations, where the nearest 
external conductor is far from the electron layer com- 
pared to the mean spacing between electrons in the layer, 
the value of | V/i| will be relatively small compared to the 
value of |V0| or jV^/ej. This is due to the fact that, 
in the absence of external screening, perturbations in 
and fi are not really independent and their characteristic 
magnitudes can be expressed via each other. Suppose 
that a non-equilibrium perturbation in the chemical po- 
tential i5/i(r) with a large length scale is created in a 
two-dimensional system. The resulting variation in the 
particle density has a magnitude Sn ~ Sfj,/{dfj,/dn). The 
magnitude of the potential variation caused by this ac- 
cumulation and depletion of electrons can be estimated 
as S4> (e^ I K)5nl^. Thus, in a compressible system, we 
find the ratio of gradients V/x/V(/) to be of the order of 
~ ~ Ts/l^ which is small in the limit of large l^. 

In an incompressible system (for instance, a macroscopi- 
cally wide strip of a system in the middle of a quantized 
Hall plateau, if there are no localized states in the energy 
gap) , the accumulation of electrons due to the perturba- 
tion in the chemical potential occurs only at the edges 
of the system. The accumulated charge creates an elec- 
tric field that slowly (as l/r) vanishes into the interior 
of the sample. The resulting ratio d^/d(j) is as small as 
l/\n{W/rs) where W is the width of the strip 

The above observation has an impact on the issue of 
the edge-versus-bulk current distribution. Since, in the 
absence of a temperature gradient, the boundary cur- 
rents arise solely due to (5/i, and the bulk currents are 
due to both dfi and 64>, the non-equilibrium current in 
a system without gates flows predominately in the bulk. 
The situation changes, however, if the temperature is not 
uniform. As we shall see in Sec. IV, the boundary frac- 
tion of the net thermocurrent is at least as significant as 
its bulk counterpart. 



III. LINEAR RESPONSE IN THE ABSENCE OF 
DISORDER 



In the previous section, the local current response to 
the electric and statistical fields was expressed in terms 
of the coefficients of the mechanical response in the bulk 
L*^*^ defined by the general expressions (p2|)-(55). Now 
we derive explicitly three of the four response coefficients 
i^'-* for a simple case of a uniform disorder-free sample. 
Since the internal magnetization currents and the bound- 
ary currents are already expressed in terms of the equi- 
librium magnetizations and their derivatives as given by 
Eqs. ( |Tl|Jl^]8l| , ^ ) , here we consider only the "transport" 
components of the bulk currents. We restrict our atten- 
tion to the case of a two-dimensional electron system, 
with magnetic field B in the z direction, perpendicular 
to the layer. 

Our main result is the expression for the number cur- 
rent density 



(83) 



which shows that, in the absence of disorder, the trans- 
port contribution to the number current density is fully 
determined by the equilibrium number and entropy den- 
sities. (In this section, the symbol e represents the two- 
dimensional antisymmetric tensor exy — —^yx — 

Equation ( |83| ) is valid locally for any interacting elec- 
tron system, provided it does not have a shear modulus, 
i.e., is a fluid, and the energy spectrum of electrons is 
quadratic. We present two different ways to obtain this 
result. In Sec. A, we derive Eq. (p3|), and hence the coef- 
ficients L^^^ and L^^^ directly, using arguments based on 
fluid dynamics. In the alternative proof given in Sec. B, 
we first derive the coefficients L^^^ and L'^^^ by study- 
ing the current response in a uniform electric field, and 
then obtain L*-^) from the Onsager symmetry. The ad- 
vantage of the second derivation is that it deals with a 
homogeneous non-equilibrium system, VT = V/i — 0. 

Both derivations that follow employ the physical no- 
tion of internal pressure, which is not quite obvious in 
the presence of a magnetic field and which we now dis- 
cuss briefly. In the presence of a uniform magnetic field 
B\\z the equation for conservation of momentum on the 
microscopic scale may be written 



m- 



d J{r) 



-p{r)VV{r) ~ -J{r) x B - V • 7T(r), 



(84) 



where V{r) is the one-body potential and 7r(r) is the in- 
ternal stress tensor at point r. As in the case of the en- 
ergy current discussed in Section II, there is not a unique 
definition of tt for a system with finite range forces. How- 
ever, a definition consistent with Eq. ( p4| ) and with the 
requirement of quasi-locality is 
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(85) 



where t^°' is defined by Eq. (29). Equation (85) may 
be checked by integrating both sides of (Q) over an in- 
finitesimal volume, to obtain the rate of change of the 
momentum inside the volume. The integral of the last 
term on the right hand side of (^J) is equal to the in- 
tegral of the stress tensor over the surface enclosing the 
volume. The first term of (|8^) then gives the change 
in momentum due to particles crossing the surface, while 
the second term gives the force exerted on particles inside 
the volume by particles outside the volume [cf . Eq. (^] . 

In applying the above equations to the present prob- 
lem, we interpret m as the band mass rather than the 
bare mass of the electron, and V{r) excludes the periodic 
potential of the ions. Thus, when there are no impurities 
present and no applied electric field, V{r) is a constant 
in the interior of the sample, and VV" arises only from 
the confining potential at the boundaries. We also as- 
sume that Uij depends only on the distance between the 
electrons, so that Fij is parallel to (r^ — rj) and the 
stress tensor is symmetric. We specialize to the case of 
a two-dimensional electron system and thus use notation 
appropriate to two dimensions. 

In thermal equilibrium, far from the boundaries, in an 
isotropic system, the stress tensor must be proportional 
to the unit tensor, so we may write 



rir)) 



(86) 



where Pint , the "internal pressure" , depends on the chem- 
ical potential and temperature. In the presence of a 
magnetic field. Pint is not equal to the pressure P which 
appeared in the thermodynamic equations (S)-(12) and 



which is equal to the force per unit length exerted by the 
boundaries on the contained electron gas. Rather, we 
have 



P, 



P - MB. 



(87) 



The difference between Pint and P arises from the Lorentz 
force exerted by B on the boundary current I — (c/ e)h x 
M. 

Equatio n (pTj ) may be obtained directly from equations 
~|) and 



( p4| ) and ( ^6D if we integrate the right-hand side of 
along a line-segment from a point ri in the interior of the 
sample, to a point r2 where p{r) = 0, passing through 
a point T on the boundary. Since dJ/dt = in equilib- 
rium, and since 7r(r) = at point r2, we find 



Pint + MP = / dr-p{r)VVir) 



The right-hand side of (88) is just the force per unit 
length exerted by the boundary at point r. The fact 
that the pressure entering the thermodynamic equations 
is indeed the same as the force per unit length exerted 
by the boundaries follows from the well-known fact that 



a magnetic field constant in time cannot produce work 
on charged particles. Hence, although there is a momen- 
tum exchange between the system and the source of the 
field via the Lorentz force, there is no energy exchange at 
B — const. Therefore, the term PdV in Eq. ( |59| ) repre- 
sents the work done by the expanding system on the ex- 
ternal confinement. To avoid confusion, we note that the 
work done by the expanding system against the Lorentz 
force actually goes to increase of the internal energy of 
the system itself pQ] . 



A. Fluid Dynamics Approach 

Consider an electron liquid in a uniform magnetic field 
and in the presence of an electric potential, chemical po- 
tential, and temperature all of which vary smoothly in 
space. All the fields and currents in the system are as- 
sumed to be either constant in time or varying at a small 
frequency as discussed in Sec. II A. We concentrate on 



a small macroscopic element of the liquid with area SA 
in the interior of the sample, which we choose to be of 
a size much less than the length scales of the fields and 
much larger than the average inter-electron distance. 

Then, setting V ^ (j) in the RHS of (|8|), and setting 
dJ/dt = 0, as is appropriate for a quasi-equilibrium sit- 
uation, we find 



V • TT + - J X B = 0, 



(89) 



where n, tt and J are averaged over the element SA. If 
the induced current is small, then the correction to tt 
arising from the current should be second order in J, 
and therefore negligible. Thus tt may be replaced by its 
equilibrium value, Pint'^a/s, evaluated for the local values 
of /i and T. To first order, Eq. (89) becomes 



J c / 1 
v = - = —e[ V(/)-H -VPint 
n eB \ n 



(90) 



Apart from the additional term resulting from the pres- 
sure gradient, the right-hand side of Eq. ( |90| ) represents 
the classical drift velocity in the crossed magnetic and 
electric fields. Substituting Eq. (|^ into Eq. ( pO| ) and 
using the relation nV/x = VP — s VP which follows from 
Eqs. (PI) and (pO|), we finally obtain 



(91) 



J = nv = Jtr - (c/e)eVM, 



where Jtv coincides with the right-hand side of Eq. 
and the second term is the internal magnetization current 
as defined by Eqs. (||) and (0). 

We note that our arguments do not apply to an elec- 
tron solid. In a solid, a non-uniform drift current will 
cause a shear deformation which will increase until the 
stress forces suppress the local drift. Equation (^) de- 
termines the force acting on an element of a liquid and 
does not include the shear stress contribution when we 
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use (86). In addition, the pinning effects which arise in 
the presence of even a weak disorder potential makes our 
consideration completely inapplicable in the case of an 
electron solid. 



B. Derivation Using Onsager Symmetry 

We now study a homogeneous system, VT = V/i = 0, 
which is driven out of equilibrium by an external electric 
field E — (l/e)V0. Since we study coefficients of the 
bulk response, we can assume that the electric field is 
uniform and the system itself is infinite. In addition to 
the original laboratory frame, we consider the system in 
the "primed" reference frame moving at a velocity 



V ^c[Ex B]/B^ = ceV0/(eB), 



(92) 



in which the applied electric field vanishes. Since the 
system is homogeneous, and disorder potential is absent, 
Galilean invariance requires that properties of the system 
in both reference frames are the exactly same. Due to the 
absence of electric field in the moving frame, the system 
with respect to this frame is in equilibrium. The number 
and energy current densities in this frame are therefore 
zero, J' = = 0. In the laboratory frame, the system 
moves as a whole at a velocity v. Hence the number cur- 
rent density defined as the average number of electrons 
passing through unit length per unit time is given by 

J = nv. (93) 

To determine the energy current density, we split the 
whole system in two parts by an imaginary straight line 
perpendicular to the drift velocity v and find out how 
much energy AE is transferred, in the laboratory frame, 
from one part of the system to another in time At. The 
energy AE has two contributions, one from the direct 
transfer of an element of the system across the line, and 
another due to the work done by one part of the system 
on the other part while moving 



AE = eAxAL + Pint AxAL, 



(94) 



where Ax — vAt is displacement of the system in time 
At, AL is the length of the line segment which we con- 
sider, and Pint is the pressure in the interior of the sam- 
ple. Substituting Eq. ( p7| ) for the internal pressure into 
(U), for the energy current density = AE/{ALAt) 
we obtain 



{e + Pi,,t)v = {e + P-MB)v. 



(95) 



One can also obtain Eq. (£^) directly from the mi- 
croscopic expression for the energy current, given by 
Eqs. (p7|)-(^. Let us write Vi = v'^ + v, where is 
the velocity of particle i in the frame moving with ve- 
locity V. In the moving frame, the system is at local 
equilibrium with an energy density eo and a stress tensor 



PintS^j. Comparing (|2|)-(|^) with (||), we see that ^) 
is correct to first order in v. 

Identifying Eqs. (^3|) and (| 9 5| ) with the expressions for 
the bulk current response (|5C|)7 (51), and substituting v 
from Eq. (p2h, wc find the response coefficients 



1(3) 



c{e + P- BM) ^ 
eB ^' 



(96) 
(97) 



As one can see from the last formula, L^^^ is an odd 
function of the magnetic field. From this fact, and from 
the symm etry relation L^^(P) — Lp^{—B) discussed in 



Sec. II El we have 

1(2) _ 1(3) _ cje + P-BM) . 

eB 



(98) 



Substituting the obtained coefficients L'^^K L'^'^^ , and 
L(3) into (|^)-(||), and using (|) and we find 

7V(i) = L(i), and 



sc , 



Then using ( |64|) we arrive at Eq. 
number current. 



(99) 

for the transport 



IV. THERMOPOWER MEASUREMENTS 

A. Thermopower Measurements and Current 
Distributions 

We now turn to discuss the thermoelectric properties 
of real samples, in which temperature gradients are main- 
tained by the coupling of the electron gas to the phonons 
of the substrate. To use the results of the previous sec- 
tions, we shall assume that the sample is homogeneous, 
and that the coupling to the substrate is sufficiently weak 
that the corresponding thermal relaxation rate is much 
slower than the microscopic relaxation rate of the elec- 
tron gas r~-^, and the response of the electron gas is well- 
described by the transport properties of the isolated elec- 
tron gas. However, we assume that the thermal coupling 
to the substrate is sufficiently strong that on a macro- 
scopic scale we may assume that the local temperature 
of the two-dimensional electron gas is equal to the local 
temperature of the substrate, and we need not impose 
the condition that V • = in the electron gas. The 
substrate is assumed to be an electrical insulator, how- 
ever, so that V • J = in the electron system. We do 
not discuss the energy current in this section. 

A convenient way in which to study the thermoelec- 
tric response of a sample is through the thermopower. A 
thermopower measurement involves the application of a 
uniform temperature gradient VT to a sample which is 
disconnected from current leads. Since there can be no 
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average electron current flow, an electrochemical poten- 
tial gradient = V(0-|-/x) develops. The thermopower 
tensor S is defined in terms of this potential gradient. 

If the diagonal matrix elements of N^^"^ and iV^^^ are 
different from zero, it can be shown that the conditions 
V • Jtr — 0, with n ■ Jti = at the sample boundaries, 
together with Eq. (^4|), require that Jtr — everywhere. 
This is the case only if 



= -SVT 

e 

everywhere in the sample, where 

5- = -(eT)-i[/V(i)]^i[7V(2)] 



(100) 



(101) 



In the special situation where TV^^^ and 7V*^^^ are pro- 
portional to the antisymmetric tensor e, as occurs for in- 
stance in the case of zero impurities, the value of ^(r) is 
not properly determined by the conservation equations, 
together with Eq. ( |64[ ) and the boundary conditions. For 
example, one may add to ^ any function f{r) which van- 
ishes at the boundaries of the sample, without affecting 
the values of V • Jtv in the interior or n • Jtr at the bound- 
ary. The average value of is still given correctly by 
Eq. (IOC), however, for any solution of the equations, and 
the value of ^ at any point of the boundary will be the 
same as if were uniform in the sample. Thus the ex- 
perimentally measured thermopower, in which the volt- 
age drop is measured between t wo p oints at the bound- 
ary, would still be given by Eq. (IOC) in this case. 

For the remainder of this section, we will focus on the 
thermopower of systems in which the disorder potential 
is weak. In the limit of vanishing disorder, one can use 
the response coefficients we have derived in the previous 
section, and the thermopower tensor takes a particularly 
simple form 



Sal3 — Saf}- 

en 



(102) 



Thus, the thermopower tensor is diagonal, with a mag- 
nitude given by the entropy per particle, s/n, divided by 
the charge per particle, — e. This result is familiar for 
non-interacting electrons ||l|,^j5| . 

Unfortunately it difficult to provide a general criterion 
for how weak the impurity scattering must be in order 
that its effects on the thermopower can be neglected and 
Eq. (102) applies. Rather, the form of such a criterion de- 
pends on the nature of the low-lying charged excitations 
of the system, which may be quite different at different 
filling fractions (compare, for example, filling fractions at 
which the system is compressible and incompressible in 
the zero- temperature limit). A necessary condition for 
Eq. ( |l02| ) to apply is that the impurity scattering is suf- 
ficiently weak that both of the tensors 7V^^^ and N^"^^ are 
almost purely off-diagonal, such that the thermopower 
tensor itself is close to diagonal. In some circumstances 
this condition may not be sufficient, as there may be cor- 
rections to the size of the diagonal thermopower. The 



form of such corrections depends on the specific experi- 
mental conditions, and requires a specific calculation of 
the effects of impurity scattering on the carriers. In the 
following, we will concentrate on thermopower of systems 
for which the impurity scattering is sufficiently weak that 
Eq. (|l02| ) applies. In the next subsection we will discuss 
the form of the corrections that can arise due to impurity 
scattering for filling fractions close to i/ = 1/2 or 3/2. 

Although, under the conditions of the experiment, no 
net current passes through the sample, in a quantizing 
magnetic field circulating non-equilibrium currents are 
induced. In the bulk, these are internal magnetization 
currents, whose continuity at the edge is provided by the 
boundary currents. As we shall now show, these currents 
can be very large, in the sense that the local current 
density in the presence of both the temperature gradient 
and the compensating electric field can be comparable to 
what it would have been in the presence of only one of 
these fields. 

Let us compare the average and the local current den- 
sities induced by a uniform temperature gradient alone, 
VT = const, with Wfj, — W(f> — 0. We consider a filling 
factor V — 1/2 which is an important and much stud- 
ied example of a compressible state. The average current 
density is equal to the transport current density (|83|) 



Jtr = ^eVT. 

eB 



(103) 



The portion of this current that fiows in the bulk is given 
by 



J bulk T 
— Jtr 



I ybulk 

*^ mag 



cs 


c ds 




'eB 


^ edB 





eVT, (104) 



where we used Eq. (48) and the thermodynamic relations 
(|)-(0)- 

In a strong magnetic field, for which all electrons 
are restricted to the lowest spin-polarized Landau level, 
one may express the entropy per unit area in the form 
s = nQSq[v, [e^ / nt) /T], where ng — eB/hc is the num- 
ber density of flux quanta, Sq is the entropy per flux 
quantum, v = n/riQ is the flUing fraction, and we have 
assumed a Coulomb force-law, for which the typical en- 
ergy scale is set by the magnetic length, £ = y/hc/eB. 
Using d\a-{e'^ / k€) / dB — 1/2, one can write 



ds 
OB 



s 
B 



n dSg T ds 



B dh 



2B dT 



(105) 



n.B 



Now, at a filling fraction of one half, particle- hole sym- 
metry requires that dSq/dv = 0. Hence Eq. (104) may 
be rewritten 



rbulk 



{u = 1/2) = - 



cT ds 



2eB dT 



eVT. 



(106) 



n.B 



We shall discuss two cases. 
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First, for an ideal non-interacting electron gas (for ex- 
ample K —^ oo), the entropy is independent of the tem- 
perature at fixed filling fraction, so the above expression 
is identically zero. The temperature gradient induces no 
number current in the bulk of the system, and all of the 
induced current flows around the edge of the sample. 

Secondly, for a system interacting by Coulomb forces, 
the entropy at = 1/2 is believed to be approximately 
linear in temperature (with a logarithmic correction at 
very low temperature due to the divergence of the ef- 
fective mass of composite fermions) pT| . The bulk cur- 
rent induced by the temperature gradient is therefore ap- 
proximately cs/(2ei?)eVr, which is one half of the total 
current induced by the temperature gradient (103); the 
remaining current flows around the edge of the sample. 

In both of the above cases the current driven by the 
temperature gradient is found to be very inhomogeneous, 
with all or half of the total current flowing on the edge 
of the sample. 

In contrast, the distribution between edge and bulk 
of the current driven by Vi^ in a thcrmopower experi- 
ment depends on the apportionment between V/i and 
V(^. The relative proportions of coming from V/x and 
'V(j> depends on the compressibility of the electron system 
and on the electrical capacitance per unit area, i.e., on 
the distance to the nearest conductivity surface. In most 
practical cases, the contribution of (/) will be much larger 
than fj,, so that V(/) ~ V^, and V// « 0. In this case, 
the compensating current driven by Vi^ will be uniformly 
spread over the sample and a strong circulatory current 
is set up by the combination of and VT. The non- 
equilibrium part of the current density induced by the 
thermopower measurement has a form shown schemati- 
cally in Figure 2. 




FIG. 2. Schematic diagram of the distribution of the addi- 
tional current induced in a thermopower measurement in the 
quantum Hall regime. No net current flows. However, a large 
fraction of the current induced by the temperature gradient 
VT is at the edges of the sample, whereas the compensating 
current induced by the electric field V<^ is spread uniformly 
over the interior. 



B. Experimental Comparison 

Finally, we will compare our conclusions concerning 
the thermopower of a sample in the limit of weak disor- 



der with recent thermopower measurements in the frac- 
tional quantum Hall regime. At high temperatures, the 
observed thermopower is dominated by the phonon-drag 
contribution resulting from the momentum exchange be- 
tween the system and the phonons in the substrate. Very 
low temperatures are required before the intrinsic ther- 
mopower caused by the diffusion and drift in the sys- 
tem itself can be observed. It is only very recently that 
this has been achieved in the fractional quantum Hall 
regime. In Ref. |^,|lO|, thermopower measurements on 
a hole gas in this regime are reported, and it is found 
that at temperatures below about 100 mK the intrinsic 
thermopower can be distinguished. The crossover from 
a dependence of the thermopower at high temper- 
atures, to a linear temperature-dependence below 100 
mK is associated with the transition from phonon-drag- 
dominated to diffusion-and-drift thermopower. At such 
low temperatures the signal-to-noise ratio is rather poor 
in thermoelectric measurements, so it is difficult to re- 
solve much structure related to the incompressible states 
at fractional filling fractions. One can, however, distin- 
guish dips in the diagonal thermopower, Sxx, at i/ = 1/3, 
2/5, 3/5, 2/3, consistent with the expectation that the 
thermopower should vanish at these filling fractions for 
which the entropy is exponentially small. 

A more interesting issue is the behavior at even- 
denominator filling fractions. It is found that Sxx ex- 
hibits a broad maximum at v — 1/2 and, less clearly, at 
I' = 3/2. As explained in Ref. M the absolute values 
of the thermopower are inconsistent with a single par- 
ticle picture, and electron-electron interactions must be 
considered. 

It has been argued that, at even denominator filling 
fractions, the appropriate description of the electron sys- 
tem is in terms of a Fermi liquid of weakly interacting 
composite fermions |2^. We will use this model in con- 
junction with equation ( ^3|) to calculate the thermopower 
of a disorder-free system. Let us first assume that elec- 
trons are maximally spin-polarized. Thus, we assume 
that at a filling fraction of ly'i = i + 1/2, the lowest i (spin- 
split) Landau levels are filled (and therefore contribute 
zero entropy), and that the remaining half- filled Landau 
level may be represented by a Fermi liquid of compos- 
ite fermions with effective mass m*. The entropy of the 
half-filled level is determined by the density of states at 
the Fermi surface, so we obtain the thermopower of the 
system to be 



^af3 - -T-— 2 J (Ja/: 

o Ti nq 



(107) 



where the particle charge q is — e for electrons and +e 
for holes. Note that the same formula applies at all half- 
integer filling fractions i^i — i + 1/2. However, one must 
be careful to use the appropriate value for the effective 
mass, which may vary for different absolute magnetic 
fields and for different filling fractions i^i. 

Corrections to Eq. (107) due to impurity scattering 
may arise if the scattering rate t^^ of the composite 
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fermions from impurities is larger than the microscopic 
equilibration rate of the disorder- free composite- 
fermion system. Since this equilibration rate becomes 
very small for composite fermions close to the Fermi sur- 
face, it may well be the case that in experiments at low- 
temperature the impurity scattering rate is suffi ciently 
large, 3> r~^, that corrections to Eq. (107) arise. 
Arguments based on a Boltzmann transport theory for 
the composite fermions suggest that, in this cas e, im- 
purity scattering will affect the thermopower ( |107| ) if the 
scattering rate is energy-dependent [Q. Specifically, 
if we consider a model of composite fermions with con- 
ventional parabolic dispersion, E (x k^, and a transport 
scattering rate that varies with energy as E~p. then the 



effect of impurity scattering is to multiply Eq. ( 107 ) by 
a factor {1 + p): 
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Q/3, 



(108) 



This has the same form as the conventional Mott for- 
mula for the thermopower of a spinless two-dimensional 
electron gas in zero magnetic field [|| . 

A calculation of the scattering of composite 

fermions in modulation-doped quantum wells suggests 
that it is only weakly energy-dependent, p ~ 0.13, 
and the corrections to (107) are small. In the follow- 
ing, we will compare the experimental observations with 
Eq. (107), bearing in mind that a prefactor (l+p) — 1.13 
may arise due to impurity scattering. 

Comparing Eq. ( |107D with the measurements of Sxx 
reported by Ying et al. |9j, we find that, at a filling frac- 
tion of h' — 1/2 and at a magnetic field B = 5.6T, an 
effective mass of m* = 1.3 ± 0.3mo is required for con- 
sistency, where mg is the free-electron mass. This value 
is a factor of 2 larger than the value m* sa 0.7mo ob- 
tained in Ref. from their own analysis of the data, 
which used the Mott formula Eq. (108) with an assumed 
value of p ~ 1. However, our value of the effective mass 
does not seem inconsistent with estimates of m* based 
on other types of transport measurements. For example, 
a value of m* — lAnio at = 1/2 is quoted in Ref. for 
a hole doped sample with a higher carrier density, such 
that f = 1/2 occurred at B = 13T. In an ideal system, 
with no Landau-level mixing, zero layer thickness, and no 
impurity corrections, the effective mass at = 1/2 would 
be expected to be proportional to \/B, which would pre- 
dict a value m* « 0.9rno for the sample used for the 
thermopower measurements, where = 1/2 occurred at 
5.6T. However, it is not at all clear that this scaling 
should apply to the actual samples. (The observed value 
of the effective mass at 13T is in any case considerably 
larger than one would expect based on numerical studies 
of finite systems where Landau-level mixing, finite thick- 
ness, and impurity effects are ignored.) 

We note that the value p ~ 1 assumed in Ref. |^ was 
obtained from calculations of impurity scattering of elec- 
trons in zero field, whereas calculations pS] for composite 



fermions suggest a much smaller value p ~ 0.13, as was 
mentioned above. 

The difference between our formula and the one used 
in Ref. |^ is much more serious at v = 3/2. In that 
case Ying et al. ^ replace the particle density n in 
Eq. (108) by the density of composite fermions, which 
is now 3 times smaller than the density of holes in the 
valence band. We, believe, however, that whether one 
uses Eq. (107) or Eq. (|108|) the quantity n should be 



the total number of carriers, including those in any filled 
Landau levels. 

Since the experimental thermopower reported in 
Ref. Q is larger at = 3/2 than at ly = 1/2 by a factor 
« 1.4, Ying et al. conclude that to*(3/2) ~ 0.5to*(1/2), 
which they view as evidence for the validity of a model 
of spin-polarized composite fermio ns a t v — 3/2. How- 
ever, we would conclude using Eq. ( |l07|) or Eq. (108) that 
m*(3/2) w 1.4m*(l/2). This is contrary to the expecta- 
tion for the ideal system that to*(3/2) should be smaller 
than m*(l/2), due to the smaller value of B. 

Unfortunately, we do not see any justification for the 
analysis used by Ying et al. a t — 3/2. We be- 



lieve that our starting formula (107) is correct in the 
limit of small impurity scattering, and that energy- 
dependent impurity scattering leads to an additional 
prefactor (l+p) that is close to unity. Moreover, the rela- 
tion TO*(3/2)/m*(l/2) = 5^^(3/2)/S'^^(l/2) holds even 
allowing such impurity scattering, provided the exponent 
p is the same at each filling fraction. (In the absence of 
any Landau level mixing this would necessarily be the 
case if the magnetic length were the same at each filling 
fraction. The change in magnetic length by a factor of 
\/3 that occurs between i/ — 1/2 and 3/2 at fixed electron 
number density is unlikely to affect the exponent p.) 

It therefore appears to us that the reported ther- 
mopower measurement at = 3/2 arc not consistent with 
a simple model based on spin-aligned composite fermions. 
The failure of this model may be due to the combined ef- 
fects of increasing degree of Landau level coupling and 
the smaller Zeeman energy expected at ;y = 3/2 as com- 
pared to V = 1/2. Alternatively, the effects of disorder 
may be quite different at these two filling fractions. Ne- 
glecting any significant effects of disorder, however, and 
viewing the diagonal thermopower as a measure of the 
entropy, one would conclude that the entropy at = 3/2 
is larger than what one would expect from a model of 
maximally spin-polarized composite fermions. It may 
be that additional entropy arises from the loss of spin- 
polarization. A number of experiments indicate that in 
typical electron doped GaAs samples, the electron sys- 
tem is not maximally spin-polarized at v = 3/2 even at 
T = ||2^-^7[| . It is not clear whether this will also occur 
for hole-doped samples, where the Zeeman energy may be 
more important due to the larger g-factor. To gain better 
understanding of the origin of the discrepancy at = 3/2, 
it would be interesting to investigate the dependence of 
the thermopower on the extent of Landau level coupling 
(e.g., by studying n-type samples, or p-type samples with 
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different densities) , and on the Zeeman energy (by tilted 
field measurements). 



V. CONCLUSIONS 



We have discussed the linear response of a homoge- 
neous, bounded interacting electron gas in quantizing 
magnetic field. We studied the number and energy cur- 
rents which arise in response to the gradients in electric 
and chemical potential and in temperature. We derived 
general expressions for the bulk and boundary currents in 
the presence of mechanical and statistical fields. In gen- 
eral, the boundary of the sample can carry a finite frac- 
tion of the total current passing through the sample. The 
local response in the bulk may be described as a sum of 
"transport" and "internal magnetization" contributions. 
Internal magnetization currents do not contribute to the 
net current, are always divergenceless, and cannot be re- 
vealed in any standard transport experiments performed 
on either homogeneous or inhomogeneous samples. They 
can be detected only in special contactless experiments 
resolving the local current distribution. We found that 
Onsager symmetry relations cannot, in general, be ap- 
plied directly to the local current densities in the bulk of 
the sample. However, they do hold locally for the trans- 
port currents, and therefore for the net currents passing 
through the sample. We derived expressions for three 
of the four response functions of an interacting system in 
the limit of weak disorder in terms of equilibrium proper- 
ties of the system. In particular, we showed that, in this 
case, the thermopower tensor is diagonal and is propor- 
tional to the entropy per particle. Recent thermopower 
measurements on a high-mobility sample show that this 
conclusion is consistent with a model of a Fermi liquid of 
spin-polarized composite- fermions at v = 1/2. However, 
for the observations to be consistent with this model at 
V = 3/2, a very large effective mass is required. An effec- 
tive mass of this size seems unlikely, and we suggest that 
the spin-polarized composite-fermion state may not be a 
good description of the system at that filling fraction. 
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